Distribution of fractal dimensions at the Anderson transition 
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We investigated numerically the distribution of participation numbers in the 3d Anderson tight- 
binding model at the localization-delocalization threshold. These numbers in one disordered system 
experience strong level- to- level fluctuations in a wide energy range. The fluctuations grow substan- 
tially with increasing size of the system. We argue that the fluctuations of the correlation dimension, 
D2 of the wave functions are the main reason for this. The distribution of these correlation dimen- 
sions at the transition is calculated. In the thermodynamic limit (L — > 00) it does not depend on 
the system size L. An interesting feature of this limiting distribution is that it vanishes exactly at 
D2 m ax = 1.83, the highest possible value of the correlation dimension at the Anderson threshold in 
this model. 
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The localization-delocalization Anderson transition 
has posed for a long time a fascinating problem. In sys- 
tems with short range interaction, purely diagonal disor- 
der and unbroken time-reversal symmetry, without spin- 
orbit interaction, it occurs for dimensions d > 2 jlj]. In 
the thermodynamic limit (i.e. infinite system size) the 
transition point separates the systems where all wave 
functions are localized from the systems where some part 
of them is extended. Exactly at the transition one finds in 
the center of the band extended wave functions. Due to 
the proximity of the unavoidable localization at one side 
of the transition they have a self-similar fractal (actually 
multifractal) structure. This is a direct consequence of 
quantum critical fluctuations at the transition point. 

This multifractality of the wave functions at the local- 
ization threshold is one of the most important features 
discovered [§J3| since the pioneering work of Anderson [[| . 
This fruitful idea has became widely recognized (see e.g. 
Refs. [||-0]) and has considerably helped our understand- 
ing of different phenomena in mesoscopic systems related 
to electron localization. For instance the well known 
log-normal distribution of the conductance in disordered 
metals || can be taken as a finger-print of multifractal 
wave functions which survive in the weakly disordered 
state (so-called pre- localized states [Q). 

Usually the multifractal (as well as fractal) structure 
of a wave function manifests itself in the size dependence 
of the participation number (PN) 
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where L is the system size and D2 < d is the correlation 
dimension of the wave function ^(r). For a localized state 
D2 = and M do not depend on L. On the other hand 
D2 = d for a delocalized wave function which extends 
uniformly over the sample. The inequality < D2 < d 



means that a multifractal wave function is delocalized 
but, in the thermodynamic limit, nevertheless occupies 
only an infinitesimal fraction of the sample. 

Due to strong level-to-level fluctuations |J it is very 
hard to verify the size dependence of M for a particu- 
lar state in a computer experiment. At the transition 
these fluctuations increase with increasing system size. 
Therefore, one has to study the size dependence of the 
averages of M or, preferably, the distribution function. 
The size dependence of the fluctuations of M is then con- 
verted to the size dependence of this distribution func- 
tion. Choosing suitable size dependent variables one can 
collapse these distributions (for different system sizes) to 
one universal curve. However, to do this systematically 
one needs to understand the origin of these fluctuations. 
In our opinion the main source of above mentioned giant 
fluctuations of PN at the transition are the fluctuations 
of the fractal dimension D2 of the wave functions. 

In a disordered system, exactly at the transition [ [l0[ , 
critical wave functions with very different degrees of de- 
localization (participation numbers) should coexist inde- 
pendent of their energy. This is a direct consequence of 
an infinitesimal neighborhood of the localized behavior 
at one side of the transition and delocalized behavior at 
another side. Accordingly each delocalized wave function 
should have a different D2. If the distribution function 
V{F>2) becomes size-independent in the thermodynamic 
limit our hypothesis will be correct. 

In the present paper we present numerical results for 
this distribution (more exactly for the distribution of log- 
arithms of the PN) and show that it is indeed universal, 
i.e. it does not depend on the system size in the ther- 
modynamic limit. This is somewhat reminiscent of the 
previous idea of Shapiro [|ll],|l2| about the existence of 
universal distributions at the Anderson transition. This 
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problem has recently been addressed analytically in a 
couple of papers, but far from the transition point. 

In Ref. [|l3| the variance of the inverse participation ra- 
tio (IPR) was calculated for disordered metallic samples 
in leading order of small parameter 1/g 2 , where g ^> 1 
is a dimensionless conductance of the system. It was 
suggested that the relative value of the IPR-fluctuations 
should be of order unity at the transition point and that 
their distribution should be universal. A similar question 
was raised recently in Ref. Q where the distribution of 
the IPR was calculated for a large but finite conductance 
g of small metallic grains. 

To investigate this problem at the transition we numer- 
ically solve the standard Anderson tight-binding model 
with diagonal disorder on a 3d simple cubic lattice with 
the Hamiltonian 



W = 5>IO(i|+£*'. 
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To model a disorder we distribute the site energies Si 
uniformly in the interval — W/2 < S{ < W/2. For the off- 
diagonal elements we took tij = 1 for nearest neighbors 
and otherwise zero. The Anderson transition is in this 
model at a critical value, W c = 16.5 (see, e.g. Ref. 
and references therein). 
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FIG. 1. Participation numbers versus energy for a 3d sys- 
tem of 14 x 14 x 14 = 2744 sites at the Anderson transition. 

Diagonalizing the matrix Hij for a cubic sample with 
N = L 3 sites and open boundary conditions we find a set 
of N orthonormal eigenvectors e s (j) and corresponding 
eigenvalues, Ej. The participation number for a state j 
is defined as usual by 



N 



(3) 



Fig. [j] shows the participation numbers Afj versus en- 
ergy Ej for a 3d system with N = 2744 sites. At the 
transition there is a strong level-to-level fluctuation of 
the PN. In the whole energy range states are found which 



are very close in energy but whose Afj differ by up to two 
orders of magnitude. This means that the energy of a 
particular state is not indicative of the spatial behavior 
of the wave function. In a next step we, therefore, cal- 
culate the distribution function of the PN in an energy 
band around zero energy (where the Anderson transition 
takes place). We take the width of this band equal 10 
i.e. we include all states with \Ej \ < 5. The distribution 
remains practically the same if we take a much narrower 
strip, \Ej\ < 1 but the computation time increases sig- 
nificantly. 
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FIG. 2. Distribution functions of participation numbers at 

the Anderson transition for different system sizes. In the 
brackets the numbers of realizations used in the averaging 
are shown 



The normalized distribution functions of the partici- 
pation numbers for different system sizes are shown in 
Figure ||. As expected the distributions are strongly size 
dependent. Increasing the system size the position of 
the maximum shifts to higher values, approximately as 
oc A/" 0,3 , and the amplitude of the maximum decreases as 
oc A?- - 46 . 

Let us suppose now that the fluctuation of the PN at 
the transition is due to fluctuations of D2 of the wave 
functions. According to Eq. (|l|), without loss of gener- 
ality, we can take this relation in the form J\f oc N D<2 ^ . 
Then if Fn(J\T) is the normalized distribution function of 
PN for an ensemble of systems with N sites and different 
disorder, we can extract the distribution function of the 
correlation dimensions, D2, in this ensemble 



V N (D 2 ) = (1/S)F N (V° 2/3 ) N D2/3 \nN. 



(4) 



Fig. [| shows that this distribution of correlation dimen- 
sions is much less sensitive to the system size than the 
one of the PN. Moreover, with increasing N it obviously 
approaches some size independent function, Voo (^2), i.e. 
in the thermodynamic limit a universal distribution func- 
tion of correlation dimensions does exist at the Anderson 
transition. For a given model of disorder it should be the 
same for different realizations. In other words we believe 
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that the distribution function Voo(D2) is a self-averaged 
quantity [|l6[ and can be obtained by analyzing one suf- 
ficiently large system. 

Fig. U shows an additional interesting feature of this 
distribution: with increasing argument the function 
Vn{D2) rapidly decays to zero. The drop gets steeper 
with increasing size N . This is seen more clearly in a 
logarithmic plot of Vn {D2). For systems with N > 512, 
Vn(D2) drops four orders of magnitude in the small in- 
terval 1.5 < D2 < 2. We think, therefore, that in the 
thermodynamic limit there is a point, £>2max where the 
function Voo(D2) approaches exactly zero. In such a case 
this would be the highest possible value of the correla- 
tion dimension D2 at the Anderson transition. In the 
following we present additional evidence in support of 
this idea. 




FIG. 3. Distribution functions of the correlation dimension 
Vn(D2) at the Anderson transition for different system sizes. 

For that purpose we calculate the size dependence of 
the average PN and of the PN of the the most extended 
state in the system, i.e. for a given realization of disor- 
der the state with maximal participation number, A/" max - 
The ensemble averages of both quantities (denoted by 
aapn and ampn, respectively) against system size N are 
plotted in Fig. |j. First, both quantities scale with TV as 
a power law, aN D2 ^ 3 . Secondly the value D2 = 1.26 for 
the average participation number (aapn) is very close to 
the position of the maximum in the distribution function 
of correlation dimensions shown in Figure ||. It is sur- 
prising that averaging over all states does not affect the 
power law behavior given by Eq. ([l]) for a single state. 
For the ensemble average of the maximum participation 
number, (A/* m ax) we find D2 — 1.83. This is the corre- 
lation dimension of the most extended state in the sys- 
tem at the transition. To answer the question whether it 
really is the maximal correlation dimension, I^max, we 
should study the fluctuations of this quantity. A maximal 
correlation dimension Z^max does exist if the fluctuation 
of this quantity goes to zero with N — > 00. 
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FIG. 4. Averaged maximum participation number (ampn) 
and averaged average participation number (aapn) at the 
transition versus system size. Lines are the best least squares 
fit with (AT) = aN D2/3 , where (AT) is "ampn" or "aapn", 
respectively. 



Fig. H shows the relative fluctuations of the maximum 
participation number versus system size. Firstly the fluc- 
tuations are rather small and decrease slowly with in- 
creasing N. According to Eq. (|l|) we can relate them to 
the fluctuations of the correlation dimension 5D2 m &x 



AAfmax/ (A/max) = (5£>2max/3) In N. 



(5) 



We conclude from the above that, in the thermodynamic 
limit (N — > 00), the fluctuation of the correlation dimen- 
sion for the most extended state SD 2max — » at least not 
slower than 1/lniV. Therefore Z^max = 1-83 is indeed 
the highest possible value of the wave function correlation 
dimension at the Anderson transition. The distribution 
function V 00(^2) should approach exactly zero at this 
value. 
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FIG. 5. Relative root mean square deviation of the max- 
imum participation number against system size. The solid 
fitting line corresponds to a power law behavior, oc 1/N ' 13 . 

There are numerous previous calculations of the corre- 
lation dimension D2 at the 3d Anderson transition. From 
the scaling with system size of the density-density corre- 
lation of the wave functions, the correlation dimension 
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was estimated to be D 2 = 1.7 ± 0.3 (for W = 16) @. 
From the size dependence of the participation number 
averaged, both over different wave functions in a small 
energy interval 0.25 around zero and over disorder (with 
a Gaussian distribution of site energies), the correlation 
dimension D 2 = 1.6 ±0.1 was estimated Jl8| . In Ref. |l9[ 
from the spectral compressibility, x °f the levels, us- 
ing x — (1 ~~ D 2 /d)/2 derived in Ref. Q the correla- 
tion dimension was estimated as D 2 = 1.4 ± 0.2. From 
the time decay of the temporal autocorrelation function, 
C(t) oc t~ D2 / 3 , the correlation dimension, D 2 = 1.5 ±0.2 
was obtained plfl . Using box-counting procedures D 2 — 
1.7 ± 0.2 H, L> 2 = 1.52 ± 0.11 H, and L> 2 = 1.68 || 
was calculated. Again using box counting and averaging 
the results over wave functions in a small energy interval, 
AE — 0.01, around zero and over five different realiza- 
tions of disorder D 2 = 1.46 was found [|l5|. Again using 
box-counting techniques in Ref. |2^], D 2 = 1.33 ± 0.02 
was obtained. 

This shows clearly the quite large uncertainty, in the 
literature, in the existing values of D 2 at the transition. 
The reason is obviously that one deals with a distribu- 
tion of this quantity. Among different characteristics of 
this distribution one can consider for example the posi- 
tion of the maximum, at about 1.3, the average value, 
D 2 = 1.26, and the correlation dimension of the most 
extended state at the transition, D 2mdiX = 1.83. 

Similarly the distribution function of the information 
dimension D\ as well as of the other generalized dimen- 
sions D q at the transition can be obtained. The results 
of this analysis will be published elsewhere. In the ther- 
modynamic limit all these distribution functions are ex- 
pected to be universal and self- aver aged quantities. Us- 
ing the Legendre transformation the multifractal spec- 
trum fj(a) for each wave function with energy Ej can 
be calculated (compare Ref [|B|). In other words critical 
wave functions at the transition should be characterized 
by a whole distribution of /(a)'s which is already a func- 
tional. In particular, the positions of the maxima, ao's 
of these functions should be distributed. 

In conclusion, we investigated numerically the distri- 
bution of the correlation fractal dimension D 2 for the 
critical wave functions at 3d Anderson transition. This 
distribution appears to be universal i.e. it no longer de- 
pends on the system size in the thermodynamic limit. 
Extrapolating these results to other fractal dimensions 
we conclude that each critical wave function should pos- 
sess its own (infinite) set of generalized fractal dimen- 
sions, Dd at the transition. And one should rather speak 
about the distribution functional of these functions. 
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